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We present microwave-frequency NbTiN resonators on silicon, systematically achieving internal quality factors 
above 1 M in the quantum regime. We use two techniques to reduce losses associated with two-level systems: 
an additional substrate surface treatment prior to NbTiN deposition to optimize the metal-substrate interface, 
and deep reactive-ion etching of the substrate to displace the substrate-vacuum interfaces away from high 
electric fields. The temperature and power dependence of resonator behavior indicate that two-level systems 
still contribute significantly to energy dissipation, suggesting that more interface optimization could further 
improve performance. 


Superconducting coplanar waveguide (CPW) mi¬ 
crowave resonators are crucial elements in photon de¬ 
tectors^, quantum-limited parametric amplifiers'’^ and 
narrow-band filters^, as well as read-out, interconnect 
and memory elements in quantum processors based on 
circuit quantum electrodynamics^. They also play a crit¬ 
ical role in hybrid devices, connecting superconducting 
circuits with micro- and nanomechanical resonators®’^ 
and solid-state spins®’®. In many quantum science and 
technology applications, resonators must operate in the 
quantum regime, requiring low temperatures to reach the 
ground state (thermal energy k^T small compared to the 
photon energy at resonance, hfr) and single-photon ex¬ 
citation levels. Under these conditions, however, inter¬ 
nal quality factors (Qis) are typically substantially lower 
than their high-temperature or high-power values. 

In the quantum regime, the dominant loss mecha¬ 
nism for high-Q superconducting resonators can be at¬ 
tributed to parasitic two-level systems (TLSs) in the di- 
electrics^®’^^. TLSs may reside in the bulk substrate^^, 
as well as in the metal-substrate, metal-vacuum and 
substrate-vacuum interfaces^®’^^^^® where electric fields 
may be large (see Ref. 19 and references therein for a 
recent review of material-related loss in superconduct¬ 
ing circuits). Interface TLSs are common by-products 
of the fabrication process, often introduced by impurities 
associated with substrate surfaces^®’®^ and etching chem¬ 
istry®^. To our knowledge, the best resonators reported 
to date^® {Qi = 1.72 M at 6 GHz) are fabricated by 
epitaxially growing aluminum on sapphire substrates fol¬ 
lowing careful surface preparation (high-temperature an¬ 
nealing in an oxygen atmosphere). For CPW resonators 
on silicon (Si) substrates, achieving Qi > 1 M in the 
quantum regime has proven challenging, with the best 
resonators reported in Ref. 23. 

In this letter, we present silicon-based, gigahertz- 
frequency CPW resonators with Q, systematically above 
1 M in the quantum regime, fabricated from niobium 
titanium nitride (NbTiN) superconducting films. This 
performance is reached by optimizing two aspects of 
the fabrication. First, the substrate surface is treated 
with hexamethyldisilazane (HMDS) immediately prior 



FIG. 1. False-colored scanning electron micrographs of CPW 
resonators and their cross sections, (a) Capacitive conpling 
of a RIE feedline to a DRIE resonator, (b) Schematic cross 
section of RIE and DRIE CPW resonators (vertical axis not 
to scale) in a NbTiN thin film (gold shading) on a Si sub¬ 
strate (grey). The treated substrate-metal interface is high¬ 
lighted in red. RIE (DRIE) etches ~ 150 nm (~ 80 yim) into 
the substrate, (c) Cross section of a DRIE resonator near 
the surface, showing the metallization (300 nm). Si underetch 
(~ 1 jam) and sidewall roughness (~ 100 nm). (d) Cross sec¬ 
tion of DRIE resonator, showing the high etch anisotropy and 
depth (~ 80 pm). 


to metal deposition to reduce losses associated with the 
metal-substrate interface. Second, we employ highly 
anisotropic deep reactive-ion etching (DRIE) of the sub¬ 
strate to displace lossy substrate-vacuum interfaces away 
from regions of high electric field. We obtain the best per¬ 
formance by combining both techniques, reaching 2 M in 
the quantum regime. 

Resonator fabrication starts with surface preparation 
of high-resistivity (p = 150 kUcm) Si (100) substrates 
for deposition of the superconducting film. We choose 
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NbTiN as the superconductor for its high-quality metal- 
vacuum interface^"^. Because we define the resonators 
via etching of the metal film rather than lift-off, we can 
apply surface treatments immediately prior to metal de¬ 
position to optimize the substrate-metal interface. The 
substrate surface quality depends on its physical condi¬ 
tion (e.g., roughness), chemical composition (e.g., pres¬ 
ence of native oxides or nitrides), and chemical cleanli¬ 
ness (e.g., absence of polymers or metals). For example, 
thin (~ 2 nm) amorphous native silicon oxide layers are 
known to host a considerable density of TLSs^^’^^. Var¬ 
ious adsorbates can in turn attach to the native oxides. 
For instance, hydroxyl groups (-OH) couple to dangling 
bonds on the substrate surface and build a layer of silanol 
groups (Si-OH). This silanol layer is also hydrophilic, and 
thus adsorbs additional water and dipolar molecules. We 
first remove lossy oxides from the surface following the 
standard practice of dipping the Si substrate in a 1:7 
buffered solution of hydrofluoric acid (HF) for 2 min, 
which both etches away surface oxides and terminates 
the Si surface with hydrogen. We then include an ad¬ 
ditional step immediately following the HF dip, placing 
the substrate on a hot plate for 2 min at 110°C while 
exposing the surface to an HMDS-nitrogen(N 2 ) atmo¬ 
sphere. HMDS is an organosilicon compound most com¬ 
monly used as a primer to improve resist adhesion, to 
create hydrophobic surfaces, and to prepare substrates 
in the growth of high-mobility graphene^^. XPS surface 
analysis of substrates immediately after these treatments 
revealed an absence of Si-0 bonds for both the HF-only 
and HF-I-HMDS cases^®. After also analyzing the dy¬ 
namics of surface re-oxidation^®, we conclude that the 
Si surface does not re-oxidize either during the HMDS 
treatment or loading of the substrate in the high-vacuum 
chamber for NbTiN deposition. Following surface prepa¬ 
ration, NbTiN is sputtered in a reactive atmosphere of 
argon (Ar) and N 2 at 3 mTorr (5 % atm -of N 2 ) and 
400 W (DC). The films are 160 or 300 nm thick^®, with a 
typical critical temperature of 15.5 K and p = 110 pHcm. 

We next aim to reduce intrinsic losses associated 
with the substrate-vacuum interface. We employ deep 
reactive-ion etching (DRIB) of the exposed Si substrate 
using the Bosch process^® to move the substrate-vacuum 
interface away from regions of high electric field (Fig. la- 
b). The Bosch process, consisting of alternating etching 
and passivation steps, is highly tunable, with variable 
etch depth, under-etch, sidewall roughness and passiva¬ 
tion^^. The DRIB resonators presented here have ~ 1 pm 
under-etch below the edges of the metal regions, side-wall 
scalloping with ^ 100 nm features (Fig. Ic), and a PTFB 
fluorocarbon layer uniformly deposited on the exposed Si 
surfaces during the etching process. 

We now describe a recipe which allows fabrication of 
both RIB and DRIB structures on the same chip. First, 
the NbTiN films are patterned using e-beam lithogra¬ 
phy and RIB to define feedlines and standard resonators. 
Next, an 8-pm-thick layer of AZ 9216 photoresist is spun, 
followed by 25 nm of sputtered tungsten and a layer of 



FIG. 2. Intrinsic quality factors of four DRIE resonators 
measured at 15 mK. Qi is extracted from measurements of 
complex-valued feedline transmission. The DRIE resonators 
(NbTiN thickness ~ 300 pm, central conductor width W = 
12 pm, and gap S = 2 pm) have fundamental frequencies 
ranging from 2.75 to 6.41 GHz. The best-fit values of A, 
Qi and Qc are used to convert the calibrated input power to 
(uph) [25]. All resonators show Qi > 1 M in the quantum 
regime, with higher QiS attainable for higher (uph) and lower 
/r. Error bars are smaller than the symbol size. 


PMMA A4. After e-beam patterning and developing the 
PMMA, a short RIB step is used to etch the tungsten, 
creating a metallic mask for a 1 h anisotropic 02 -plasma 
ashing of the AZ resist. We then deep-etch the Si sub¬ 
strate, alternating between etching using an inductively- 
coupled plasma (ICP) of SBe gas (20 mTorr, 2200 W, 7 s) 
and passivation with an ICP of C 4 F 8 (2 mTorr, 1200 W, 
2 s). The thick AZ layer allows etching to a depth of 
80 pm (Pig. Id) and 500 nm minimum feature size. 

For the devices in this study, we employ quarter-wave 
CPW resonators with center conductor width fixed to 
IF = 12 pm and separation S between the center con¬ 
ductor and ground planes set to give a transmission¬ 
line impedance close to 50 H. For RIB resonators, we 
set S' = 5 pm, taking into account er = 11.9 for the 
silicon substrate and a kinetic inductance fraction of 
~ 16%. For DRIB resonators, we decreased S to 2 pm 
to compensate for the reduced effective dielectric con¬ 
stant. Prom resonator frequencies, we estimate a DRIB 
resonator impedance of ~ 62 H. The resonators are ca- 
pacitively coupled to a common feedline with a target 
coupling quality factor Qc between 700 K and 1 M. Af¬ 
ter dicing and final cleaning in NMP, the devices are 
mounted in a copper PCB with Apiezon N grease, wire- 
bonded and anchored to the mixing chamber plate of a 
dilution refrigerator routinely employed for characteriz¬ 
ing resonators (and qubits) in the quantum regime^®. 

We quantify resonator performance by measuring their 
power and temperature-dependent resonance lineshapes. 
Specifically, we measure the complex-valued feedline 
transmission iS '21 near the fundamental frequencies A of 
each coupled resonator, and fit the data^®’®® to extract 
/r and Qi- Here, we report measurements for our best 
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resonators, fabricated using both surface treatment and 
DRIE techniques (see supplementary material^^ for other 
devices). We show Q, as a function of excitation power 
for four HMDS-DRIE resonators with fundamental fre¬ 
quencies /r = 2.75 — 6.41 GHz (Fig. 2). All resonators 
show Qi > 1 M in the quantum regime (average intrares¬ 
onator photon number (riph) ~ 1, T ~ 15 mK). However, 
we observe further increases in Qi (up to 4 times at 2.75 
GHz) as the resonator is populated with more photons^®. 
This suggests that the low-power limits for the resonators 
are still significantly influenced by TLSs. We also observe 
a systematic increase in Qi with decreasing /r, explana¬ 
tions for which include the decrease in TLS polarization 
at lower frequencies and the decrease in radiation losses 
at longer wavelengths. 

To study the effects of TLS polarization, we acquire 
the power dependent Qi of the 2.75 GHz resonator for 
T = 15-400 mK (Fig. 3a). For all excitation powers, 
Qi increases monotonically with temperature over this 
range, and higher temperatures make Qi less power de¬ 
pendent, with the factor of 4 change at 15 mK dropping 
to a factor of only 1.5 at around 400 mK. We therefore 
conclude that TLSs still provide significant loss in the 
quantum regime, and attribute the increase of Qi with 
temperature to thermal depolarization of TLSs. 

A more detailed temperature dependence is acquired 
by measuring S '21 around 2.75 GHz at fixed {ripQ ~ 10® 
while the device is continuously cooled from 900 to 
15 mK. Fitting S '21 as before, we extract Qi and /r. 
Initially, Qi increases during cooling due to decreasing 
quasiparticle-induced loss^®, before reaching a maximum 
value of ~ 8.6 M at T ~ 0.4 K (Fig. 3b). Below 0.4 K, 
Qi decreases again as increasing TLS polarization leads 
to greater absorption of resonator photons. 

Resonator loss from TLSs predominantly depends on 
their density of states at /r and the applied microwave 
power^®’®®. The continuing power dependence of Qi at 
single-photon levels (Fig. 2) shows that even at these low 
powers, resonant TLSs are still being depolarized by the 
incident microwaves, suggesting a significant loss contri¬ 
bution from TLSs with low saturation powers or, equiv¬ 
alently, long lifetimes. Because Qi does not reach its 
low-power baseline, this limits the usefulness of the Qi 
measurements in quantifying the TLS-related loss. An al¬ 
ternative way to probe the effect of TLSs is to study the 
temperature dependence of /r. Unlike loss, the disper¬ 
sive shift of the resonator also depends on non-resonant 
TLSs that are not saturated by the applied microwaves®^. 
Their influence on /j- can therefore be measured even at 
high (uph). Figure 3c shows the temperature dependence 
of the resonator frequency shift Afr{T) = /r(r) — /r(0) 
extracted from the same data as Fig. 3b. The local max¬ 
imum in A/r at r ~ 500 mK again reflects the loss- 
dependent pull of the resonator, but unlike Qi, A/r is 
not monotonic below this temperature. 

The low-temperature behavior of A/^ is well described 
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FIG. 3. Temperature dependence of 2.75 GHz DRIE res¬ 
onator. (a) Power dependence of Qi measured for T = 15- 
400 mK. At high powers, Qi is similar for all T, but temper¬ 
ature dependence increases for lower (riph). (b,c) Resonator 
probed at fixed {tipQ ~ 10® (black arrow in (a)) while device 
is cooled from 900 to 15 mK (plotted against mixing chamber 
temperature). Both Qi (b) and A/ (c) show non-monotonic 
temperature dependence. We attribute the initial increase 
in Qi and A/ (as temperature decreases) to a decrease in 
quasiparticle-related loss. At T < 400 mK, Qi and A/ begin 
to decrease as TLS-related losses increase. Data htting to a 
simple model®® for A/ (solid red line in (c); Eq. 1) shows 
that, in this regime, the behaviour is dominated by the T- 
dependence of TLS polarization. Modifying the related TLS 
model for Qi to account for high excitation power (solid red 
line in (b); Eq. 2) provides estimates for a TLS polarization 
fraction at (riph) ~ 10® of ~ 4 %, and a background (non- 
TLS) loss of ~ 1/(9 M) which is consistent with high-power 
QiS from (a). 


by a two-parameter TLS model®® 
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Here, fk is the complex digamma function, and /r(0) 
is the resonator frequency at T = 0, and k = FS is 
a combined loss parameter, where S is the loss tangent 
(at T = 0 and (riph) = 0) and F the filling factor of 
the TLS host material. Fitting the TLS model with k 
and /r(0) as free parameters (red curve in Fig. 3c) shows 
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good agreement with the data for T < 400 mK and ac¬ 
curately captures the distinctive local minimum in A/ at 
T ~ 60 mK. The best-fit value k = 6.35(4) • 10“^ com¬ 
pares well with previous values found for Nb resonators 
on sapphire^®. This value can also be used to predict the 
low-power, low-temperature limit of TLS-related loss^°, 
giving Qi « 1/fc ~ 1.6 M. The ~ 20 % disparity with 
the measured single-photon Qi (Fig. 2) fits with previous 
observations^®. It arises partly because A/r is sensitive 
to TLSs over a broad frequency bandwidth, giving an es¬ 
timate of average loss tangent 5 which may differ from 
the local loss tangent inside the resonator bandwidth®®, 
and partly because we do not reach a baseline value of 
Qi even at the single-photon level. 

We also fit an heuristic model for the low-temperature 
dependence of Qi modified from Ref. 30: 


1 


kes tanh 


( hfrim 

V S/cbT J 


+ 


1 


Qother 


( 2 ) 


Here, we introduce two ad hoc fitting parameters, a 
modified loss parameter fceff to account for the reduced 
TLS polarization at high excitation power, and an ex¬ 
tra, fixed loss term l/Qother which accounts for the non- 
TLS-related loss which becomes significant as tempera¬ 
ture increases. This fit also agrees well with the data at 
T < 0.4 K (red curve in Fig. 3b), with best-fit values 
of fceff = 2.64(4) • 10-® and Qother = 9.02(2) • 10®. The 
ratio kes/kfit ~ 0.04 provides an estimate of the TLS po¬ 
larization fraction still remaining at (nph) ^ 10®, while 
Qother is Consistent with the high-power QjS measured in 
Fig. 3a. 

In conclusion, we have demonstrated few-GHz CPW 
resonators on silicon with intrinsic quality factors 
exceeding one million in the quantum regime, achieved 
through surface treatment and substrate etching. Our 
results suggest that these resonators exhibit a signihcant 
reduction in TLS-related losses, but the temperature and 
power dependences indicate that TLSs are still an impor¬ 
tant source of energy loss. This offers the opportunity for 
further improvements by interface optimization, perhaps 
pushing quantum-regime resonator performance towards 
the high-power limit with QjS approaching 10 M at lower 
frequencies. Since we have only partially explored the 
available parameter ranges in this study, future efforts 
will include optimizing the DRIE recipe and studying the 
effects of etch depth, side-wall roughness and under-etch 
on resonator losses. Finally, these surface-treatment and 
DRIE techniques can be directly applied to the fabri¬ 
cation of superconducting qubits in integrated circuits, 
offering a new path towards longer qubit coherence times. 
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SUPPLEMENTAL MATERIAL 

This supplement provides experimental details and data sets supporting claims made in the main text. First, we 
present design and fabrication details for three device generations that we have measured: Device C (reported in 
the main text) and devices from Sets A and B (reported below). We then describe several technical details of the 
measurement system: the measurement setup and device shielding, the fitting of the complex-valued transmission 
spectra, and the procedure for estimating (riph). Next, we present measurement results for devices from Sets A and B. 
Finally, we discuss the results of XPS analysis of the pre-deposition substrate surface following different treatments. 


DEVICES 



FIG. SI. a) Optical image of HMDS-treated device from Set B (Set B used the same layout as Set A). The feedline and the 
lower 6 resonators are etched by RIE while the top 6 resonators are patterned by DRIE. Freqnencies for the RIE and DRIE 
resonators are centered around 7 and 4 GHz, respectively, b) Optical image of Device C (HMDS+DRIE). The feedline is 
patterned using RIE, bnt all resonators are fabricated using DRIE with frequencies spanning 2 to 11 GHz. 


Devices in Sets A and B had identical layouts (Fig. Sla). Each device contained 6 RIE resonators with frequencies 
~ 7 GHz (spaced by ~ 100 MHz; see results in Figs. S4 and S6) and 6 DRIE resonators with frequencies ~ 4 GHz 
(also spaced by ^ 100 MHz; see results in Figs. S5 and S7). In each set, we fabricated one device using the standard 
HF-only substrate treatment (green dashed curves in Figs. S4-S7) and one device using the HF-I-HMDS substrate 
treatment (solid red curves in Figs. S4-S7). In Sets A and B, the DRIE recipe is slightly different from the recipe 
described in the main text (for Device C). For these sets, the plasma power used in the DRIE etching step was set 
at 1800 W. At the end of the deep etch process, the samples underwent an additional etch with an SFg inductively 
coupled plasma (20 mTorr, 1800 W, 60 s with no alternating passivation step). This removed the PTFE passivation 
and created a slightly isotropic etch profile. The fabrication was otherwise almost identical for Sets A and B, with 
the NbTiN film ~ 160 nm thick and the DRIE etch ^ 80 /im deep. The only difference was that the devices in Set B 
were fabricated with the inclusion of a short (20 sec) additional fluorocarbon (PTFE) side-wall passivation step in 
the DRIE recipe after the isotropic etch, while no such final step was included in the devices from Set A. 

Device G (recipe and data in main text) only included DRIE resonators fabricated using the combined HF-l-HMDS 
surface treatment. We used a thicker (^ 300 nm) NbTiN film, and kept the DRIE etch depth the same (~ 80 /rm). 
The device had 10 resonators with frequencies between 2 and 11 GHz. The plasma power in the etching step in the 
DRIE recipe was increased to 2200 W to reduce the feature size of the side-wall scalloping. 
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MEASUREMENT SETUP 



FIG. S2. Measurement setup used for resonator characterization in the ®He/'‘He dilution refrigerator (Leiden Cryogenics CS81): 
microwave schematic, magnetic shields (two Cryophy, one aluminum), IR Liters (home-made Eccosorb filters), and IR adsorbers 
(inner surface of the inner Cu radiation shield, coated with Stycast 2850 and silicon carbide powder). 

In an effort to minimize resonator losses induced by non-equilibrium quasiparticles®^ and magnetic vortex displace¬ 
ment®^, we mount the samples inside several layers of shielding (Fig. S2). Specifically, we anchor the PCB-mounted 
sample directly to a copper cold finger connected to the mixing chamber of the dilution refrigerator. The sample is 
then enclosed in a copper can, the inner surface of which is coated in a mixture of Stycast 2850 and silicon carbide 
granules with diameters between 15 and 1000 /im. This can is enclosed in a second can made from aluminum. This is 
finally enclosed in two layers of cryogenic magnetic shielding (1-mm-thick Cryophy, Magnetic Shields Ltd). Coaxial 
cables enter the sample chamber through tubulations in the lid of the outer layer of magnetic shielding to reduce 
the impact of open holes on the shielding effectiveness. Extra radiation shielding is provided by home-made in-line 
Eccosorb filters in the input coaxial line mounted outside the magnetic shields at the mixing chamber stage. 


Data fitting 


The transmission through the full measurement chain near one resonance is shown in Fig. S3. We fit the S '21 curves 
following the method presented in Ref. S3. Equation 11 in Ref. S3 describes the general form of a microwave resonator 
transmission measurement with a number of non-ideal elements, including both inductive and capacitive coupling, 
and impedance mismatch in the feedline: 


S'21 


/ r'‘ 

V 1 + 2zQi 



(S3) 


Here, /r is the resonance frequency, A is the transmission amplitude away from resonance, Qi is the loaded quality 
factor of the resonator, related to Qi and Qc by l/Q\ = l/Qc + 1/Qi, and Qe = IQel exp(—fd) is a complex-valued 
quality factor related to Qc via 1/Qc = Re(l/Qe). The imaginary part of 1/Qe gives rise to an asymmetry in the 
resonator lineshape. We further modify Eq. S3 as follows: 
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FIG. S3. Typical measured response of a CPW resonator around its resonance frequency, and fit. Data taken from DRIE 
resonator with /r = 3.68 GHz at (riph) ~ 10 photons (Fig. 3 of main text). Panel (a) shows the trajectory of S 21 in the complex 
plane with a fit (red line) of Eq. S4. (b) Normalized |52r| calculated from the data and fit in (a). 


Here, a allows a linear variation in the overall transmission chain in the narrow frequency range around any given 
resonance, and and (j)^ account for the propagation delay to and from the sample. An example best fit of Eq. S4 
to data is also shown in Fig. S3. 


Calibration of the average number of photons in the resonator 


To estimate (uph), we calibrated the signal input line and the home-made eccosorb filter at low temperature in 
order to estimate the power Pin at the feedline input. To calculate (uph) from Pin, we note that the reflected and 
transmitted powers (ignoring feedline mismatch) are respectively given by: 

-Prefl =Hin|5ll|", 

-ftrans = -Pin|5'21 P? 


where S 21 = Qc/iQc + Qi) and S'!! = 521 — 1 = —Qi/iQc + Qi) are the usual scattering parameters for a resonator- 
shunted feedline probed near resonance®^. The power absorbed in the resonator is therefore: 


Pahs = Pil 


Preii T^trans — 


2Q? 

QcQi 


^abs can also be written using the definition of internal loss rate K[: 


-fabs — ('^ph) ^ I 

- (nph) . 

Combining these equations then gives the relation between (uph) and Pin: 

/ 

“ hf Qe “■ 

Due to the difficulty in accounting for all reflections in the full measurement line, these values can be interpreted as 
an estimated upper limit to (uph). 











RESULTS FROM SETS A AND B AT 20 rtiK 



FIG. S4. Comparison between HF and HMDS surface 
treatment of RIE resonators from Set A. 



FIG. S6. Comparison between HF and HMDS surface 
treatment of RIE resonators from Set B. 



FIG. S5. Comparison between HF and HMDS surface 
treatment of DRIE resonators from Set A. 



FIG. S7. Comparison between HF and HMDS surface 
treatment of DRIE resonators from Set B. 
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XPS DATA 

We performed XPS surface analysis of Si substrates in order to estimate surface oxygen uptake during the HMDS 
treatment following the HF dip, as well as during the approximately two-minute interval between finishing the surface 
preparation and loading the sample into the high-vacuum chamber (Phase = 2 x 10“® Torr) for NbTiN deposition. 

Figure S8a shows that the standard HF surface treatment already removes all signatures of SiO from the substrate 
surface. No discernible change is observed when substrates are placed in a mixed nitrogen-HMDS atmosphere, created 
by HMDS bubbling, following the HF dip. Figure S8b shows the XPS spectrum after subsequent exposures to air. 
While measurement directly after fast loading does not reveal a Si-0 peak, subsequent exposures of the sample to 
atmospheric conditions show the onset of one, indicating oxygen adsorption by the substrate surface. 




FIG. S8. XPS analysis after fast loading of substrates (~ 2 min after surface processing), (a) XPS spectrum of the Si-0 peak 
at 103.5 ± 0.3 eV for: an untreated Si surface with native oxides (green), a Si surface treated by HF dip only (black), and a Si 
surface treated by HF dip and HMDS bubbling (red) as described in the main text, (b) Time evolution of the XPS Si-0 peak 
for a Si surface treated by HF dip and HMDS bubbling after subsequent exposures to atmospheric conditions (times indicated 
in the legend). (We thank V. de Rooij and E. van Veldhoven for carrying out the XPS measurements.) 
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